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We report on a search for the doubly Cabibbo suppressed decay D'^tt", based on a 

data sample of 657 x 10^ BB pairs collected at the T(4S) resonance with the Belle detector at the 
KEKB asymmetric-energy e^e~ collider. We find no significant signal and set an upper limit of 
S(-B^ D*^7r°) < 3.6 X 10^^ at the 90% confidence level. This limit can be used to constrain the 
ratio between suppressed and favored B — > D*n decay amplitudes, r < 0.051, at the 90% confidence 
level. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh, 14.40.Nd 



In the Standard Model, CP violation arises from 
a complex phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix [l|, Q. Precise measure- 
ments of CKM matrix parameters are therefore of funda- 
mental importance for the description of the weak inter- 
action of quarks and the investigation for the new sources 
of CP violation. Measurements of the time-dependent 
decay rates of B'^{B^) — > D*^Tr^ provide a theoretically 
clean method for extracting sin(20i -I- (ps) where 
and 03 arc the interior angles of the CKM triangle [J|. 
The CP violation parameters are given by Q 



5± = 



2(-l)^rsin(2(/)i +(/)3 ±<5) 
l + r2 



(1) 



where r is the ratio of the amplitudes of the doubly 
Cabibbo suppressed decay (DCSD), Z)*+7r^ to the 

Cabibbo favored decay (CFD), B" D*-ti+ (Fig. 1), L 
denotes the angular momentum of the final state, and 5 
is the strong phase difference between DCSD and CFD. 
It is difhcult to determine r from i?" decays because the 
DCSD amplitude is small compared to the contribution 
from mixing followed by CFD, B^ ^ D*+7r-. 

Using available branching fraction measurements, r 
can be expressed as 



tan 9, 



(2) 



where Oc is the Cabibbo angle, and the decay constants 
Jd* and /d* are available from lattice QCD calculations. 
However, the assumption of SU(3) symmetry and addi- 
tional W^-exchange contributions result in an uncertainty 
of about 30% on r. In order to avoid this uncertainty, 
one can instead use the isospin relation. 



/ rgo 2B{B+ D*+nO) 
tb+ B{B^ D*-TT+) ' 



(3) 



where tb+/tb'> = 1.071 ± 0.009 and B{B^ -> D*-tt+) = 
(2.76 ± 0.21) X 10-3 0. We naively estimate B{B+ 



£)*+^o-j _ g g ^ 7^ taking into account the r factor 
of 0.02 calculated from Eq. ^ 0]. The previous search 
gives an upper limit of B{B+ D*+Tr°) < 1.7 x 10"^ at 
the 90% confidence level [ll. 
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FIG. 1: Feynman tree diagrams for (a) CFD B° D*-tv^ 
with the CKM coupling KlKd, and (b) DCSD B+^°'> - 
£)*+jj.o{-) ^j^]^ ^jjg coupling V*f,Vcd- 



In this paper, we report on a search for B^ _D*+7r*' 
based on a data sample of 605 fb^^ corresponding to 
(657 ±9) X 10^ BB events, collected with the Belle detec- 
tor at the KEKB asymmetric-energy e+e" collider [l3| 
operating at the T(45') resonance. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector, a 50- 
layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel- like ar- 
rangement of time-of- flight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a magnetic field of 1.5 T. An iron flux- 
return located outside of the coil is instrumented to de- 
tect Kl mesons and to identify muons. 

To search for B^ Z3*"'"7r°, we reconstruct D*^ can- 
didates by pairing a low momentum charged pion {tt'^^^) 
and a D", which is reconstructed through its decays to 



7r+, and A'Stt+tt . Inclusion 



K-TT+, K-Tr+TT°, K 
of charge conjugate modes is implied throughout this pa- 
per. 
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For charged kaon and pion candidates except pions 
from Kg's, we require tracks to have a distance of clos- 
est approach to the interaction point within 5 cm along 
the z-axis (anti-parallel to the positron beam direction) 
and within 2 cm in a plane perpendicular to the z-axis. 
Particle identification (PID) is based on the likelihoods 
TZ{K/tt) = C.kI{C,k + C^-n)-, where Lk (-C^) is the likeli- 
hood of kaons (pions) derived from the TOF, ACC, and 
dE/dx measurements in the CDC. The PID selections, 
which are '1Z{K/tt) > 0.3 (< 0.3) for kaons (pions) are 
applied to all charged particles except pions from Kg's. 
The PID efficiencies are 94% (91%) for kaons (pions), 
while the probability of misidentifying a pion as a kaon 
(a kaon as a pion) is 12% (6%). 

Neutral pions are formed from photon pairs with an 
invariant mass between 0.118 GeV/c^ and 0.150 GeV/c^, 
corresponding to ±3 standard deviations (a). The pho- 
ton momenta are then recalculated with a tt" mass con- 
straint. We require the tt^ momentum to be greater than 
0.2 GeV/c in the center-of-mass system (c.m.s.), and the 
photon energy to be greater than 0.1 GeV in the labora- 
tory frame. 

Kg candidates are reconstructed from pion pairs of 
oppositely-charged tracks with an invariant mass be- 
tween 0.485 GeV/c^ and 0.510 GeV/c^, corresponding to 
±5ct. Each candidate must have a displaced vertex with 
a flight direction consistent with that of a Kg meson 
originating from the interaction point. Mass- and vertex- 
constrained fits are applied to obtain the 4-momcnta of 
Kg candidates. 

For selection, the invariant mass of the daughter 
particles is required to be within 3a from the nominal 
mass, where cr (~ 5 MeV/c^) depends on the decay mode. 
Z?*^ candidates are required to have a mass difference 
AM = Mdt, — M]j within 3cr from the nominal mass 
difference, where a 0.5MeV/c^) depends on the decay 
mode. Mass- and vertex-constrained fits are applied to 
Z?" and candidates. 

Wc reconstruct a candidate from a and 
a 7r° candidate. We identify B decays based on re- 
quirements on the energy difference AE = Ei — 
Zi'bcam and the beam-energy constrained mass Mbc = 

is the beam energy. 



v/^bca^m- IE^pTP, where E^o 
and Pi and Ei are the momenta and energies of the 
daughters of the reconstructed B meson candidate, all 
in the c.m.s. We select candidates in a fit region de- 
fined as |A£;| < 0.25 GeV and 5.20GeV/c2 < Mbc < 
5.29GcV/c^. The signal region is defined as \AE\ < 
0.1 GeV and 5.27GeV/c2 < Afbc < 5.29GeV/c2. 

To suppress the background from continuum (e+e^ 
qq, q ~ u,d,s,c) events, we calculate modified Fox- 
Wolfram moments [lH and combine them into a Fisher 
discriminant. We calculate a probability density func- 
tion (PDF) for this discriminant and multiply it by PDFs 
for cos 6b, Az, and cos 9h, where Ob is the polar angle 
between the B direction and the beam direction in the 



c.m.s., Az is the displacement along the beam axis be- 
tween the signal B vertex and that of the other B, and 
9h is the angle between the tt^^^ direction and the oppo- 
site of the B momentum in the D*'^ frame. The PDFs 
for signal, generic B events and continuum are obtained 
from GEANT3-based ^ Monte Carlo (MC) simulation. 
These PDFs are combined into a signal (background) 
likelihood variable i3sig(bkg) ; we then impose requirements 
on the likelihood ratio TZ = Csig/ (^Csig-l-^bkg). Additional 
background suppression is achieved through the use of a 
-B-flavor tagging algorithm [ist , which provides a discrete 
variable indicating the fiavor of the tagging B meson and 
a quality parameter ^tag, with continuous values ranging 
from for no fiavor information to unity for unambigu- 
ous fiavor assignment. The backgrounds from continuum 
and generic B events are reduced by applying a selection 
requirement on TZ for events in each rtag region that max- 
imizes the value of Nsig/ Ngig + N^kg, where A'sig and 
A^bkg denote the expected signal and background yields 
in the signal region, based on MC simulation. This re- 
quirement eliminates 99% (94%) of the background from 
continuum {B decays) in the signal region, while retain- 
ing 35% of the signal. 

The fraction of events with more than one candidate 
is 3%. We select the best D*~^Tr° candidate based on the 
value of x?ot = xL(r>o)+XAM + xL(^o), where eachx^ is 
defined as the squared ratio of the deviation of the mea- 
sured parameter from the expected signal value and the 
corresponding resolution. The reconstruction efficiency is 
determined to be 0.56%, using the fitting procedure de- 
scribed below for the signal MC samples. The branching 
fractions of D*~^ and are included in the efficiency 

After the selection criteria are applied, the dominant 
background sources in the fit region are the continuum 
events and B'^ D*~^p~, while other B decays such 
as B^ and B^ L'*°7r*' have smaller contribu- 

tions. To obtain the signal yield, we perform an unbinned 
two-dimensional (2D) extended-maximum-likelihood fit 
to the AE-Mhc distributions in the fit region. The likeli- 
hood function consists of the following components: sig- 
nal, continuum background {qq), B'^ D*^ , and 
other B decays. 

The likelihood function for the signal is defined sepa- 
rately for each of the four decay modes and unified 
using the available branching fractions of the D° sub- 
decays 0, while those for qq and backgrounds from B 
decays are defined as the sum of four decay modes. 
Each AE and Mbc shape for the signal is modeled by the 
sum of a Gaussian and a bifurcated Gaussian with means 
and widths fixed to the values obtained from MC simu- 
lation. The AE and Mbc PDFs for qq are modeled by a 
linear function and an ARGUS function 1J| , respectively. 



The backgrounds from B ' — + D*^ p~ and other B decays 
are modeled by the superposition of Gaussian distribu- 
tions constructed from unbinned MC events, where the 
width of each Gaussian represents the smoothing param- 
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eter for the event [l5|. The D*^p background 

forms a large peak in the region AE < — O.lGcV and 
5.27GeV/c2 < Mbc < 5.29GeV/c2. The size and shape 
of the D*^p~ component strongly depend on the 

fraction of the longitudinal helicity component (|i?o|); we 
use |ffo| = 0.941 from Ref. 

The following parameters are allowed to vary: qq PDF 
parameters and yields of signal, qq and — > D*^ 
components. The yield of other B decays is fixed to the 
branching fractions in Ref. Q . 

Figure [5] shows the results of the fit to the data in the 
fit region. The projections of the fitted B signal in AE 
(Mbc) in the Mbc i^E) signal region are shown. We 
obtain 4.5^3 4 B~^ _D*+7r° signal candidates in the 
signal region (statistical error only). The significance is 
1.4(7, defined by ^-2 ln(£o/>Cmax) where £,„ax (Co) is 
the likelihood value at the maximum (with the signal 
fixed to zero). The likelihood function is convolved with 
an asymmetric Gaussian distribution that represents the 
systematic error. 
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FIG. 2: Projections of the unbinned two-dimensional like- 
lihood fit to data in the region |A_B| < 0.25 GeV and 
5.20GeV/c^ < Mbc < 5.29GeV/c^ (a) AE distribution for 
5.27GeV/c^ < Mbc < 5.29GeV/c^ with a magnified view 
of \AE\ < 0.07 GeV in the inset, (b) Mbc distribution for 
|A_E| < 0.1 GeV. The points with error bars represent the 
data, while the curves represent the various components from 
the fit: signal (thick solid line), continuum (dash-dotted line), 
D*^ p~ decay (dotted line), other B decays (dashed 
line), and the sum of all components (thin solid line). 



The systematic error components proportional to the 
signal yield are determined as follows. We estimate 
the systematic error from the TZ requirement by apply- 
ing the TZ requirement to data and MC events using a 
B~ D'^p^ control sample. The systematic error on the 
AA'I requirement is estimated by applying the AM re- 
quirement to W D*+7r- data and B+ D*+7r° MC 
samples. The systematic error on the secondary branch- 
ing fraction is calculated from errors given in Ref. 
The systematic error due to the charged-track recon- 
struction efficiency is estimated to be 1.0% (1.6%) per 



charged kaon (pion) using partially reconstructed _D*+ 
events. The systematic error due to TZ{K/tt) selection has 
a relative uncertainty of 0.8% (1.4%) per charged kaon 
(pion), determined from D*+ D°Tr+, D° K-tt+ 
decays. The tt" reconstruction is verified by comparing 
the ratio of £>° — > K^n^ and /v ^tt+tt" yields with 

the MC expectation; an uncertainty of 3.0% per particle 
is assigned. The Kg reconstruction is verified by com- 
paring the ratio of Z)+ ^s''^^ ^'^^ /v^7r+7r+ 
yields with the MC expectation; an uncertainty of 4.9% 
is assigned. The systematic error due to the signal MC 
statistics is 0.5% and the error due to the uncertainty 
in the total number of BB pairs is 1.4%. The system- 
atic error components proportional to the signal yield are 
summarized in Table HI 

The systematic errors on the yield extraction are esti- 
mated as follows. We estimate the uncertainty of \Ho \ of 
5° D*^ p^ by varying | Jfo| by zblcr, where the error 
of \Hq\ is taken from Ref. [lg|. Possible AE shifts be- 
tween data and MC simulation for the B^ — > D*^p~ 
background are evaluated by measuring the AE' shift 
of the B^ D*'^p~ background component using a 
B'^ D*'^tt'^ control sample. To obtain the systematic 
error on the background fraction of other B decays, we 
vary the normalizations of the individual sources by zblcr, 
where the values are taken from Ref. The normal- 
ization of other background components are varied by 
±50%. The systematic error due to the uncertainty in 
the shape of the B background PDF is determined by 
varying the Gaussian smoothing width by factors of two 
and one half from its nominal value. Uncertainties from 
the two-dimensional correlation in the signal and the qq 
components are estimated by applying 2D background 
PDFs to the signal and the qq shapes. The effect of a 
possible bias in the fitting procedure is estimated by a 
toy MC study. The systematic errors on the yield ex- 
traction in the signal region are summarized in Table ITIl 



TABLE I: Systematic errors for B{B^ 
tional to the signal yield. 



D* 



prop or- 



Source 


Systematic error (%) 




±a 


TZ requirement 


3.0 


AM requirement 


3.3 


Secondary branching fractions 


3.3 


Track finding efiiciency 


5.1 


Particle identification 


4.4 


tt" reconstruction 


4.1 


Kg reconstruction 


0.3 


MC statistics 


0.5 


Number of BB pairs 


1.4 


Quadratic sum 


9.8 



We then obtain the branching fraction of i?+ 
tohc B{B+ D*+Tr°) = [1.2;j-](stat)t°j(syst)]xl0- 
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TABLE II: Systematic errors for B{B+ - 
the yield extraction in the signal region. 



i3*+7r°), related to 



Source 



Systematic error 
(number of events) 





+ 0" 


— a 


\Ho\ of ^ D'+p- 


0.7 


-1.9 


AE shift oiW D*+p~ 


0.0 


-0.6 


Fraction of backgrounds 


0.8 


-0.4 


Gaussian width of 


0.5 


-2.0 


background PDF 






2D correlation for qq 


0.0 


-1.3 


and B+ D*+t7° 






Fit bias 


0.0 


-0.5 


Quadratic sum 


1.2 


-3.2 



The likelihood distribution (£), which is convolved 
with the systematic error, is used to obtain the upper 
limit on the branching fraction. We calculate the 90% 
confidence level (C.L.) upper limit (UL) using the rela- 
CdB = 0.9 to be 



tion J^^ CdB/ 

B{B+ 



D 



< 3.6 X 10^ 



(4) 



The obtained upper limit is consistent with the naive 
estimate, 5.9 x 10~^ discussed above. This result can be 
used to obtain an upper limit on the ratio of magnitudes 
of DCSD and CFD in D*tt decay. 



r < 0.051 (90% C.L.). 



(5) 



To summarize, a search for the doubly Cabibbo sup- 
pressed decay D*^Tr'^ in a data sample of 605 fb~^ 
yields an upper limit of B{B+ D*+Tr°) < 3.6 x 10"^ 
at the 90% confidence level. This limit can be used 
to constrain the ratio between suppressed and favored 
B D*Tr decay amplitudes, r < 0.051, at the 90% con- 
fidence level. 
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